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ABSTRACT 

We present the results of Rossi X-ray Timing Explorer (RXTE) and Swift 
monitoring observations of the magnetar IE 1547.0—5408 following the pulsar's 
radiative outbursts in 2008 October and 2009 January. We report on a study of 
the evolution of the timing properties and the pulsed flux from 2008 October 4 
through 2009 December 26. In our timing study, a phase- coherent analysis shows 
that for the first 29 days following the 2008 outburst, there was a very fast in- 
crease in the magnitude of the rotational frequency derivative z>, such that i) was 
a factor of ~ 60 larger than that reported in data from 2007. This z> magnitude 
increase occurred in concert with the decay of the pulsed flux following the start 
of the 2008 event. Following the 2009 outburst, for the first 23 days, i) was con- 
sistent with zero, and v had returned to close to its 2007 value. In contrast to 
the 2008 event, the 2009 outburst showed a major increase in persistent flux, rel- 
atively little change in the pulsed flux, and sudden significant spectral hardening 
~15 days after the outburst. We show that, excluding the month following each 
of the outbursts, and because of the noise and the sparsity in the data, multiple 
plausible timing solutions fit the pulsar's frequency behavior. We note similari- 
ties in the behavior of IE 1547.0—5408 following the 2008 outburst to that seen 
in the AXP IE 1048.1-5937 following its 2001-2002 outburst and discuss this in 
terms of the magnetar model. 
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Introduction 



Anomalous X-ray Pulsars (AXPs) are young, isolated pulsars generally with a large 
inferred magnetic fields (> 10 14 G). They are detected across the electromagnetic spectrum 
from the radio band (in 3 cases) to the hard X-ray regime. Like a closely related class 
of pulsars, the Soft Gamma Repeaters (SGRs), AXPs exhibit a wide range of variability, 
including but not limited to flux changes, spectral variability, X-ray bursts, X-ray flares, pulse 
profile changes, timing glitches, and rapid changes in the rotational frequency derivative. 
A XPs and SGRs are genera l ly believed to be forms of magnet ars. The magne tar model 



( IThompson fc Duncan! Il995l ; iThompson fc Duncan! Il996l ; iThompson et al.l 120021 ) identifies 



the power source of these objects to be the decay of their strong magnetic fields. Most of 
the known magnetars have been observed to have entered an active phase at least once in 



reviews, see 



Kaspi 


(2007 


), 


Mereghetti 


(2008 


), or 



The X-ray source IE 1547.0—54083 was discovered in 1980 wi th the Einstein satellit e 
during a search for X-ray counterparts of u nidentified X-ray sources (|Lamb fc Markertlll98ll ). 
It was first proposed to be a magnetar by iGelfand fc Gaenslerl (120071 ) who showed that the 
spectrum was well described by a blackbody plus power-law model, typical of AXPs, and that 
the source is located at the center of a radio shell which is possibly a previously unidentified 
supernova remnant. 



Since its discovery, the source has exhibited many X-ray flux variations (see lBernardini et al 

201 ll . for a review). First, a comparison of the X-ray flux of IE 1547.0 —5408 between 1980 



and 2 006 observed by Einstein, ASCA, XMM, and Chandra done by IGelfand &: Gaensler 
( 120071 ) revealed that the source's absorbed 0.5—10 keV X-ray flux decreased during this pe- 
ri od by a factor of ~ 7. Then, analysis of Swift data taken in 2007 June and presented 

^3 times higher than the 



m 



Camilo et al.l (120071 ) showed that the source's X-ray flux was 
highest historic level and ~16 times higher than the lowest level. 



Camilo et al.l ( 120071 ) observed IE 1547.0—5408 in the radio band and detected pulsations 
at a period of ~ 2.069 s, the smallest yet seen in a magnetar. The detected pulse profile 
was much wider than those of most rotation-powered long-period pulsars, and it was highly 
variable. Camilo et al. also estimated P to be 2xl0 -11 , which, together with the value of 
P, imply a magnetar strength field (assuming B = 3.2 x 10 19 GV PP) of 2.2 x 10 14 G), and 
furthermore that this is the magnetar with the largest @ |z>| yet known. 



1 k\so known as SGR 1550-5418. 

2 See the online magnetar catalog at www. physics. mcgill.ca/^pulsar/magnetar/main. html. 
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Following the discovery of radio pulsations, lHalpern et al.l ( 120081 ) studied the X-ray 
emission of IE 1547.0-5408 from 2007 June to 2007 October with Swift and XMM. They 
showed that the X-ray flux faded by about ~ 50% from June to August. It subsequently 
stabilized. This suggests that an X-ray outburst occurred prior to 2007 June. Halpern et 
al. also detected faint X-ray pulsations for the first time with XMM. Moreover, during 
this period of time, a radio timing ana lysis showed that the pulsar's v had systematically 
decreased by ~25% JCamilo et al.ll2008h . 



In 2008 October, IE 1547.0—5408 entered a new active phase. The Swift BAT in- 
strument detected several soft (< 100 keV) bursts from the location of IE 1547.0—5408 
(jlsrael et al.ll2010l ). An analysis of X-ray data from a variety of telescopes following the 
burst showed that the source was >50 times brighter than the lowest reported level and tha t 
the spectrum was significantly harder (jlsrael et al.ll2010l ; lKaneko et al.ll2010t iNg et al.ll201ll ). 
The data further showed X-ray pulsations and a frequency derivative more negative -and 
changing - compared with the value reported in the radio data in 2007. 



The source entered a second, even more active phase in 2009 January (IMereghetti et al. 



(Tieneo et al. 


2010 


)■ 


Enoto et al. 


(2010a) 



reported the dis covery of a hard X-ray tail just after the outburst using Suzaku, a result 
corroborated by iBernardini et al.l (120111 ). A deta iled analysis of the bursts detected with 



Swift after 2009 January outburst is presented by lScholz &: Kaspil (120111 ). in addition to an 
analysis of the Swift persistent and pulsed flux data taken immediately following the initial 
trigger, showing that the source's flux rose for 6 hrs following the trigger. 

A comparison of the 2008 and 2009 events showed interesting differences. The former 
showed a large pulsed flux increase whereas the latter, thou gh having an eno rmous persistent 
flux increase, showed only a modest pulsed flux increase. INg et al.l (120111 ) pointed out, in 
part using data reported on in this work, that the spectral variations in Chandra observations 
of the two epochs were similar, in contrast to very diffe rent i> behaviors (see ^3]). They also 
noted a significant pulsed fraction/flux anti-correlation. IBernardini et al.l ( 1201 ll ) reported on 
data from INTEGRAL, XMM-Newton, and the Chandra X-ray Observatory, finding similar 
v, flux and spectral results. They summarize the source's flux behavior since 1980 and argue 
it shows three well-defined flux 'states.' 

Here we present a detailed analysis of RXTE monitoring observations of IE 1547.0—5408 
following its 2008 October and 2009 January outbursts, supplemented by multiple Swift ob- 
servations obtained during the same time period. We report the results of an in-depth 
analysis of the timing behavior, pulsed flux changes, and pulse profile variations. Our ob- 
servations are described in Section 2. Our timing, pulse morphology, and pulsed flux and 
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hardness ratio analyses are presented, respectively, in Sections 3, 4, and 5. 



2. Observations and Pre- Analysis 



2.1. RXTE Observations 



The results presented here were obtained using the proportional counter array (PCA) on 
board RXTE. The PCA consists of an array of five collimated xenon/methane multi-anode 
proportional counter units (PCUs) operating in the 2—60 keV ran ge, with a total effec tive 
area of approximately 6500 cm 2 and a field of view of ~1° FWHM (IJahoda et al.l 119961 ) . 



For this paper, we used 114 RXTE observations. Of these, 49 observations were taken 
between 2008 October 4 and 2009 January 14, following the 2008 October outburst, and 65 
observations were taken between 2009 January 22 and 2009 December 26, following the 2009 
January outburst. The length of the observations varied between 0.6 ks and 17 ks, but most 
were between 1 ks and 6 ks. The epochs of the observations are shown in Figure [TJ Notice 
the high density of observations in the first few days following each outburst. 

Table 1 contains a list of the 49 observations obtained after the 2008 October outburst, 
showing the observation date, length, number of PCUs on, the type of analysis that each 
observation w as used for, and the number of bursts found us ing the burst s earch algorithm 
introduced in iGavriil et al.l (120021 ) and discussed further in iGavriil et al.l (120041 ) . Table 2 
contains a similar list of the 65 observations subsequent to the 2009 January outburst. 

For all observations, data were collected using the GoodXenon mode. This data mode 
records photon arrival times with l-/zs resolution and bins photon energies into one of 256 
channels. To maximize the signal-to-noise ratio, we analysed only the events from the top 
Xenon layer of each PCU. 

For each observation, we adjusted photo n arrival times at each epoc h to the solar system 
barycenter, using the source position given in lGelfand fc Gaenslerl (120071 ) (aj2ooo=15:50:54.11 
<5j2ooo =— 54:18:23.7). Then, for each observation, we extracted a time series binned at a res- 
olution of 1/32 s, in the 2—6.5 keV range (to maximize the signal to noise ratio), which 
utilized all PCUs, to be used in the timing and pulse profile analysis. Also for each observa- 
tion, we extracted time series in various energy ranges, binned at the same time resolution, 
but excluding PCUs and 1 because of the loss of their propane layers, to be used for the 
pulsed flux, pulse profile, and hardness ratio analyses. For each of the generated time series, 
7-s long intervals of data were removed around large bursts. We verified that 7-s was suffi- 
ciently long. A large burst was defined to be any with count rate greater than or equal to 
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50 counts s 1 PCU 1 at the peak. 

2.2. Swift Observations 

In this paper, Swift observations were used uniquely for timing purposes. We used all 
16 observations taken between 2008 October 03 and 2009 January 12, 14 of which were taken 
in the month following the onset of the 2008 outburst, and the remaining 2 taken in 2009 
January, see Figure [TJ We also used 8 Swift observations taken in the 23 days following the 
onset of the 2009 outburst, and all observations taken between 2009 February 13 and 2009 
December 26 that were long enough to yield a clear pulse profile. Table 3 contains a list 
of the Swift observations used in this work, all of which were in Windowed Timing (WT) 
mode. 

For each observation, we extracted photons from a 40-pixel-long strip centered on the 
source, and we adjusted the photon arrival times of each epoch to the solar system barycenter. 
Then, for each observation, we extracted a time series binned at a resolution of 1/32 s in the 
2—7 keV range, to match the RXTE band. 

3. Phase-Coherent Timing 

Time series for all 95 RXTE observations having a sufficiently high signal-to-noise ratio 
were divided in half and treated as two separate time series. Time series for all 32 Swift 
observations were not divided. Each time series was epoch-folded using an ephemeris deter- 
mined iteratively by maintaining phase coherence as we describe below. When an ephemeris 
was not available, we folded the time series using a frequency obtained from a periodogram. 
Resulting pulse profiles, with 64 phase bins, were cross-correlated in the Fourier domain 
with a high signal-to-noise-ratio template created by adding phase-aligned profiles. The 
cross-correlation returned an average pulse time-of-arrival (TOA) for each observation, cor- 
responding to a fixed pulse phase. 

The pulse phase at any time t can usually be expressed as a Taylor expansion, 

<t>{t) = M*o) + Mt ~ to) + ^o(* - to? + -M* - t ) 3 + • • ., (1) 

2 6 

where v = 1/P is the pulse frequency, v = du/dt, etc., and subscript "0" denotes a parameter 
evaluated at the reference epoch t = to- 
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To obtain ephemerides for the time periods following each outburst, we fitted the TO As 
to the above polynomial using the pulsar timing software package TEMPcH. Since the spin- 
down of this source was unstable, and since consecutive observations were not always suf- 
ficiently close to each other, phase coherence could only be unambiguously maintained for 
a period of 29 days following the first outburst, and 23 days following the second outburst. 
The two timing solutions found are represented by solid red lines in the first two panels of 



Figure El al ong with the radio timing parameters of the same source from 2007 (jCamilo et al. 



20071 . 120081 ) . The spin parameters for these two timing solutions are shown in Table 4. The 



2008 and 200 9 ephemerides ar e cons istent within 2a with those reported in llsrael et al. 



( I2010l ) and in iBernardini et al.l ( 120111 ) , although they extend for a few days longer than in 



the previous analyses. 

We made various attempts at finding the pulsar frequency outside of these time periods. 
First, for each instance where closely spaced RXTE observations were available, a frequency 
measurement was found by doing TEMPO fits through the closely spaced TOAs. These 
frequency measurements are shown as dark blue points in the first panel of Figure [2j 

Then, all 29 observations longer than 5.5 ks were broken into 2 or 3 fragments, TOAs 
were extracted for each fragment, and a frequency measurement was found for each long 
observation by fitting a linear trend to the phases with TEMPO. These frequency measure- 
ments are shown as light blue points in the first panel of Figure |2j 

We also made several attempts at finding phase-coherent timing solutions outside the 
one month following each outburst. We broke the list of all TOAs (one or two TOAs per 
observation) into 69 overlapping sets of TOAs. For each set of TOAs, we searched for a 
phase- coherent timing solution by running a TEMPO fit for every plausible combination of 
v and v. The results of the search of one such set of TOAs are presented in Figure [3J For 
each guess combination of v and z>, TEMPO returned a best-fit v and z>, and timing residuals 
from which a reduced x 2 was extracted and represented as a shade in Figure [3j with the 
darker shade corresponding to a smaller reduced x 2 ■ Note the appearance of several "islands" 
on the Figure. The center of each island corresponds to the best-fit v and z> returned by 
TEMPO. 

The fact that many islands appear in the Figure indicates that any timing solution 
found to cover the time period spanned by the set of TOAs is not unique. Instead, the 
various islands correspond to several plausible timing solutions, differing from each other by 
a small number of pulse phases. 



See http : / / www. atnf . csiro . au/research/pulsar / tempo 
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Once a shaded z/ and z> map, like that shown in Figure |3l was obtained for each overlap- 
ping set of TOAs, we collapsed each map, in the horizontal direction, into a single shaded 
column. The shade at each height of the column corresponds to the darkest shade in the 
corresponding row of the original v and z> map before the collapse. We then plotted all the 
shaded columns in Figure HI along with the various phase-connected timing solutions and 
along with the individual frequency measurements. 

For the columns near the outbursts, the number of dark spots in each shaded column 
was small, indicating a small number of possible phase-connected timing solutions. However, 
the fact that the number of dark regions in each shaded column increased as we got further 
from the outbursts clearly indicates that timing solutions found outside of the time periods 
covering the month following each outburst are not unique. Some of the longest-running 
non-unique solutions are shown as dotted red lines in Figures [2] and HI In the 2008 event, 
the loss of phase- coherence is almost certainly due to enhanced timing noise in the pulsar as 
the frequency of observations had not significantly dropped when phase coherence was lost; 
in the 2009 event, the time interval between observations was somewhat larger when phase 
coherence was lost, although the coincidence of the loss and the decay of the persistent flux 
suggests the source behavior plays an important role. 



4. Pulse Profile Analysis 

First, for each RXTE observation, we folded the data in the 2-6.5 keV band using the 
best-fit frequency found in the timing analysis (Section 3). We then cross-correlated the 
resulting profile with a standard template in order to obtain phase-aligned profiles. 

We then divided the period following the first outburst (3 October to 22 January) into 
intervals with similar net exposures, shown with the letters A through K at the bottom of 
Figure [5j We did the same for the period following the second outburst (22 January to 
26 December 26); the intervals are shown with the letters L through P at the bottom of 
Figure El 

For each time interval, we summed the aligned profiles, subtracted the DC component, 
and scaled the resulting profile so that the value of the highest bin is unity and the lowest 
point is zero. The results are presented in Figure with the time intervals marked in the 
top left corner of each profile. The different profile qualities are due to the changes in the 
pulsed flux of the pulsar. Nevertheless, clear changes in pulse morphology are apparent. For 
example, the profile peak is clearly narrower immediately following the second 2008 pulsed 
flux enhancement, and shows the full width, but is more flat-topped following the 2009 
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outburst. The re s ults s hown here are generally consistent with those of iNg et al.l ( 120111 ) and 
Bernardini et al.l ( 1201 ll ). 



Note however that in no case are any pulse profile changes so large that the profile peak, 
the assumed fiducial point, cannot be unambiguously identified. This is in contrast to other 
examples of AXP outbursts i n which multiple pe aks appeared, rendering unambiguous pulse 
numbering problematic (e.g. IWoods et al.l 120111 ) . This gives us confidence that the timing 



analysis reported in the previous Section, in particular the difficulties phase-connecting the 
timing data, are not significantly affected by pulse profile changes. Indeed the epochs in 
which the profile changes are largest (namely immediately following the outbursts) are those 
which we have successfully phase-connected. 

To study the evolution of the pulse profile with energy, we repeated the above analysis 
for energy bands 2—4 keV, 4—10 keV, and 2—10 keV. Figure M shows the evolution of the 
pulsed flux with energy for time intervals M and N indicated at the bottom of Figure [5j 
Note that the location of the peak of the profile seems to be energy-dependent, with the 
peak more to the left at higher energies, and more to the right at lower energies. This 
implies that spectral changes might cause slight shifts in the location of the pulse peak in 
the energy band used for the timing analysis (2—6.5 keV), but again, not large enough to 
cause a loss in pulse counting. Also note that pulse profiles at higher energies are not shown 
because of the very poor signal-to-noise ratio at these energies. 



5. Pulsed Flux Time Series and Hardness Ratios 

To obtain the pulsed flux for IE 1547.0—5408, for each observation, we folded the data 
and extracted aligned pulse profiles in several energy bands. For each folded profile, we 
calculated the RMS pulsed flux, 



Frms 



\ 



2Y l ((a t * + k 2 ) -«< + <>)), (2) 



k=l 



where a k is the k^ even Fourier component defined as a k = jjY^H=iPi cos (^Tvki/N), a ak 2 is 
the variance of a k , 6 fc is the odd k^ 1 Fourier component defined as b k = jjYliLiPi s ^ n {2nki/N), 
a^ 2 is the variance of b k , i refers to the phase bin, N is the total number of phase bins, Pi 
is the count rate in the i phase bin of the pulse profile, and n is the maximum number of 
Fourier harmonics used; here n—4. 

In the top panel of Figured we show the RMS pulsed flux in 2—10 keV, 2—4 keV, and 
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in 4—10 keV. Note the 2—10 keV time series was previously shown by Ng et al. (2011) and 
Scholz & Kaspi (2011). The panel shows a slow decrease in the pulsed flux after the first 
outburst interrupted by a spike 9 days after the outburst. The spike was followed by a slow 
decrease. The source behaves differently after the second outburst, when the persistent flux 
increases dramatically (Figure (2H) while the pulsed flux does not vary significantly, implying 
a large de crease in t he pu ls ed fraction follow i ng th e seco nd outburst. Th is is co nsistent with 
results of Eg et all (boilh . iBernardini et all (boilh . and lScholz fc Kaspil (boilh . 



In the bottom panel of Figure [7J we show the hardness ratio as a function of time, 
computed from ratio of the pulsed fluxes in the energy ranges 4—20 keV to 2—4 keV. The 
Figure shows a clear discontinuity in the pulsed hardness ratio ~ 15 days after the outburst, 
in the time period bracketed by the last two dotted lines, which, notably, is subsequent to 
the bulk of the pulsed flux decay, but while the persistent flux was still relatively high. We 
also computed the hardness ratio in various other energy bands, but the discontinuity was 
most significa nt for the energy range shown. Unfortunately Suzaku (jEnoto et al. Il2010ai ) and 
INTEGRAL (IBernardini et al.l 120111 ) hard X-ray observations of the source both occurred 
prior to the hardness change seen with RXTE. 



6. Discussion 



The behavior of IE 1547.0—5408 over the past several years has been extremely complex 
indeed mirroring the complexity of the behaviors of AXPs in outburst in general. In a 
collection of detailed analyses of various X-ray and gamma-ray data sets (llsrael et al.l 12010 



Kaneko et al.ll20ld : iNg et al.l 120 111 : IScholz & Kaspi 



2011; Bernardini et al. 



201 If ), including 



the present RXTE study reported on herein, a picture emerges of dramatic and varied 
multiwavelength behavior, including at least two very different outbursts from this intriguing 
source, albeit bearing some common properties. In the latter category is t he generic X-ray 
hardness /flux correlation s een now commonly in magnetar outbursts (see iRea &: Esposito 
201 ll ; IScholz &: Kaspil 120 111 for reviews), X-ray pulse profile changes, and ubiquitous X-ray 
bursting near outburst epochs, and significant changes in spin-down rate at and surrounding 
outburst epochs. We focus on this latter point here. 

Magnetar rotational evolution during and near radiative outbursts has s hown a wide 



range of behaviours, including large spin-up glitches with unusual reco veries (e.g. lKaspi et al. 



20031 ). likely spin-down events (e.g. iPalmerl l2002t [Woods et al.ll2002f). initial spi n- up events 



paired with large over-recovery, yielding n et spin-down (|Liyingstone et al.ll2010f). and large 
v variations on week-to-month time scales ( Gavriil fc Kaspi 20041 ; Dib et al. 2009 ). We note 
that throughout, no AXP has ever been observed to spin up for an extended period of 
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any duration; only changes in magnitude in the spin- down rate punctuated by occasional, 
sudden, unresolved spin-up events have yet been observed. This is in st ark contrast to 
behavior generally seen in accreting systems (e.g. IChakrabarty et al.lll997l ). Nevertheless, 



the diversity in magnetar rotational evolutions near the epochs of radiative outbursts has 
been perplexing. 

However, we note that the behavior of IE 1547.0—5408 following its 2008 outburst is 
reminiscent of unusual timin g evolution seen thus far uniquely in AXP IE 1048.1—5937 
following its 2002 flux 'flare' JGavriil k Kaspil hooi bib et al.l 120091 1 . We compare the two 
next. 

In the case of IE 1048.1—5937, the pulsar's pulsed flux had been relatively stable in sev- 
eral years of RXTE monitoring, when a small (factor of ~2) pulsed flux in crease was observed 
near the end of 2001 (see F ig. 8e), accompa nied by short SGR-like bursts (jGavriil et al.ll2002l ) 
and pulse profile changes ( iDib et al.l 120091 ). The pulsed flux relaxed back to near quiescence 
in ~3 weeks, but then there was a second, larger pulsed flux increase, this time by a factor 
of ~3 over the span of a few weeks. This larger pulsed flux 'flare' subsequently decayed on 
a time scale of several months. This latter decay occurred just prior to an increase in the 
magnitude of z>, by over an order of magnitude (Fig. 8c). During both flares, the pulsar 
suffered significant pulse profile changes, the largest of which were at the flare peak, with 
the pulse relaxing back on approximately the same time scale as the pulsed flux dropped. 
Subsequently, there was a ~2-yr period of large variations in z>, which while always remain- 
ing negati ye, varied by as much as an order of magnitude on time scales of a few weeks (see 
Fig. 15 in IDib et al.l 120091 ; Fig. 8c). Phase connection in this extremely noisy interval was 
difficult, but enabled by extremely frequent RXTE observations - approximately three per 
7-10-day interval. 

The 2008 outburst of IE 1547.0—5408 showed similar behavior (Fig. 8). The pulsed 
flux was clearly enhanced in our RXTE observations, which were triggered by SGR-like 
bursts. This enhancement then started to relax, but the decay was interrupted by a second 
pulsed flux increase ~2 weeks later, this time by a factor of ~4, similar to the increase in 
the pulsed flux in the second 2002 flare of IE 1048.1—5937 (compare Fig. 8e and f). The 
pulse profile changed, particularly near the peak of the second 'flare.' Moreover, z>, as for 
IE 1048.1—5937, dropped by a factor of ~4, but roughly in concert with the pulsed flux 
decrease, before phase-coherence was lost (Fig. 8d). It is plausible that large i> variations 
ensued, but we cannot unambiguously confirm this as the spacing of subsequent RXTE and 
Swift observations was not as dense as for IE 1048.1—5937. Clearly though, v recovered 
to its 2007 value by the time of the 2009 outburst, so another dramatic change must have 
occurred. 
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Thus, the behavior of IE 1547.0—5408 during its 2008 outburst is similar to that of 
IE 1048.1—5937 in 2001-2002, as suggested in Figure 8: in both cases there was an initial, 
smaller flare accompanied by bursting, and a much larger flare a few weeks later. In both 
cases, the flares were accompanied by pulse profile changes in which the largest changes 
occurred near the peaks of the flared). Then, importantly, a pulsed flux decay occurred near 
or simultaneous with an increase in the magnitude of i> by a factor of ~4, followed by a 
very noisy period, in which regular phase- coherent timing was impossible. We note these 
similarities exist in spite of the fact that IE 1048.1—5937 is a persistent, bright AXP, while 
IE 1547.0—5408 has the properties of a transient AXP. 

Note however that this behavior was not replicated in the 200 9 IE 1547.0—5408 ou tburst: 
in that case z> remained close to that reported in the radio data by lCamilo et al.l (120071 ). Also, 
in the 2009 outburst, the pulsed flux was barely affected, while the persistent flux increased 
by an order of magnitude (Fig. 2). Still, phase-coherence became again impossible some 
23 days after the 2009 outburst. We note low-level but significant pulsed flux variations 
subsequently, which could be relevant. 

The similarity between the 2008 outburst of IE 1547.0-5408 and the 2001-2002 outburst 
of IE 1048.1—5937 is potentially interesting. In the magnetar model, large flux enhancements 
are easy to envision given the great stresses placed on the crust due to internal magnetic field 
decay. Moreover, field-line distortions in the magnetosphere are expected as the crust shifts 
and the magnetosp here twists, hence pulse profile changes are also generically expected. 
Beloborodovl ( 20091 ) argues that magnetospheric twists will impact the pulsar spin-down in 
two possible ways. One is that a strong twist results in an enh ancement of the poloid al field 



lines, which leads to a stronger spin-down torque on the star ([Thompson et al. 



20021). This 



enhanced v is then expected to decay back to the standard value. The second way is that 
once the twist angle if) reaches a maximum value, it "boils over," with energy that would 
have been stored in the magnetosphere had if) continued to increase being carried away by 
an intermittent magnetic outflow, which can carry away significant angular momentum. In 
this picture, any delay between such intermittent magnetic outflows and the outburst onset 
is due to the time it takes if) to reach its maxim al value, though the increase in z> is expected 
to commence much earlier (iBeloborodovl 120091 ) . 



While the torque variations picture put forth by IBeloborodovl (120091 ) may provide a qual- 
itative explanation of the behavior we have seen in the 2001-2002 outburst of IE 1048.1—5937 
and the 2008 outburst of IE 1547.0—5408, it remains to be seen if the model can explain 



4 We note that the 2007 flare of IE 1048.1—5937 also sho wed a double-pe aked pulsed flux variation 
accompanied by significant pulse profile changes; see Fig. 15 in bib et all (l2009h . 
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quantitatively the large v variations observed unambiguously for IE 1048.1—5937, namely 
order-of-magnitude magnitude changes. Indeed the magnitude of the initial torque increase 
for strong (ip > 1) twists was not estimated by Beloborodov (2009), as it requires a full 
nonlinear calculation. An initial estimate for moderate twists suggests the fractional change 
in effective dipole field due to poloidal inflation varies as ip 2 , which could be checked if many 
similar outbursts are observed. 

However, in the case of IE 1547.0—5408, it is difficult to understand why the magnitude 
of v remained small and unchanged in the 23 days following the 2009 event trigger, while the 
persistent X-ray luminosity decayed from an order-of-magnitude increase. That there was 
very little change in the pulsed flux seems relevant, suggesting any poloidal field enhancement 
has to be accompanied by a pulsed flux increase. This seems reasonable if the pulsed emission 
emerges from the polar region, a hypothesis that may be testable via X-ray polarimetry. 
Perhaps then the 2009 event involved a much greater fraction of the neutron star, and, at 
least initially, impacted the poloidal region minimally. 



As discussed by lGavriil fc Kaspil ( 120041 ). disk models appear to have difficulty explaining 
the changes in pulsed flux and torques we have observ ed. This is because, at least for the 
propeller torque prescription of IChatterjee et al.l (120001 ). the X-ray luminosity is expected to 



vary as z> 7 / 3 , which is not consistent with the data. Moreover, the disparate behaviors in 
2008 and 2009 a lso pose a major cha llenge to any disk model, though are also not obviously 
explained in the iBeloborodovl (120091 ) picture. 



The long-term timing behavior of the source is summarized in Figures 2(a) and (b). 
Clearly the enhancement in the magnitude of z> is an interesting event in the pulsar's long- 
term spin-down history. The early 2007-2008 radio data (shown in green in Fig. 2a) sh o w tha t 
a quadratic polynomial was required to describe the spin-down ( jCamilo et al.l 120071 120081 ). 
with v ~60 times smaller than that seen in 2008. If the radio ephemeris is extrapolated 
to the start of the 2008 outburst, it over-predicts the frequency by approximately ~30 /is, 
consistent with there having been other episodes of enhanced spin-down. If these were small 
in number, we might have expected accompanying radiative events as in 2008, but none are 
reported. This supports the 2008 enhanced spin-down being particularly notable. 



It was noted in iDib et al.l (120091 ) that every observed AXP flare or outburst thus far were 
accompanied by a timing event, possibly caused by some unpinning of fluid vorti ces, which 
i n tur n was caused by crustal movement due to a twist propagating outward. IDib et al. 
(120091 ) also noted that the converse is not true: many AXP glitches appear to be relatively 
silent. This leads to the question of whether the onsets of the 2008 and 2009 outbursts 
of IE 1547.0—5408 were accompanied by glitches. For the 2008 outburst, this question is 
impossible to answer given the lack of observations prior to the outburst. The question is also 
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difficult to answer for the 2009 outburst given the sparsity of data and loss of phase coherence 
prior to the outburst. To provide an answer, we found a phase-coherent solution covering 
only the three observations just prior to the onset of the 2009 outburst (2 RXTE observations 
and one Swift observation, ephemeris covering MJD 54839—54845, see the small red lines 
just prior to the 2009 oubutst in Panels 1 and 2 of Figure |2]). This returned the following 



parameters for reference epoch 54845.0: z/=0. 4826070(3) s 



v- 



-5.6(1. 3)xl0 



~ 12 s~ 2 . When 



this solution is extrapolated forward in time to the date of the first post-onset observation, 
the frequency obtained is 9.6(5) xl0~ 6 s _1 larger than the known frequency at that epoch. 
This change in frequency is in a direction opposite to that which a spin-up glitch would 
cause. This does not rule out the possibili ty of a glitch, because a glitch followed by an 
over-recovery has been reported previously ([Livingstone et al.ll2010f ). 



The hardness ratio increase shown in Figure 7 is interesting. lEnoto et al. I (j2010ai ) and 
Bernardini et al.l ( 120111 ) both reported on hard X-ray observations made within 10 days after 
the 2009 outburst trigger. They showed that in addition to a soft blackbody, a single power 
law of index T ~ 1.5 could describe the data up to 200 keV. The hardness ratio as measured 
by RXTE increased significantly on MJD 54853, ~5 days after the INTEGR AL observation, 
at an epoch when the 1-10 keV emission showed no remarkable change (IScholz fc Kaspi 
20111 ). This suggests that after MJD 54853, either a single power law no longer adequately 
described the data, or the power-law index hardened. This may indicate decoupling between 
the soft and hard emission processes; hard X-ray monitoring after a future outburst could 
check this. Decoupling between the soft and hard X-ray emission mechanisms would be 
interesting, as it would argue against a causal relationship in an observed spectral co rrelation 
between these emission components (IKaspi fc Boydstun!l2010l ; lEnoto et al. Il2010bl ). 



In summary, we have reported on detailed RXTE and Swift monitoring of the AXP IE 
1547.0—5408 from 2008 October through 2009 December. We have presented pulsed flux 
time series in multiple energy bands, as well as as a detailed analysis of the source's timing 
behavior, with phase-coherence being maintained only for two ~3 wk intervals following the 

2008 and 2009 outbursts. We showed that outside this time period several plausible timing 
solutions fit the pulsar's frequency behavior in spite of relatively dense monitoring. We find 
a remarkable decrease in the magnitude of z> in concert with the falling of the pulsed flux 
following the 2008 outburst, and note a similarity in the events of the 2001-2002 outburst 
of IE 1048.1—5937. However this interesting behavior is not fully replicated following the 

2009 outburst of IE 1547.0—5408. We note a significant hardening of the source ~15 days 
after the 2009 outburst below 20 keV, at an epoch when the 1-10 keV X-ray flux shows no 
obvious change. 



The diversity of timing behaviors of AXPs and magnetars in general near the epochs 
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of radiative outbursts con tinues to be chal l engin g to understand, although t he magnet ar 
frame work first outlined by lThompson et all (120021 ) and subsequently refined by lBeloborodov 
( 120091 ) seems promising. 



We thank the RXTE and Swift staff for their efficiency in scheduling the above- described 
observations. VMK holds the Lome Trottier Chair in Astrophysics and Cosmology, and a 
Canada Research Chair, a Killam Research Fellowship, and acknowledges additional support 
from an NSERC Discovery Grant, from FQRNT via le Centre de Recherche Astrophysique 
du Quebec and the Canadian Institute for Advanced Research. 
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Table 1. RXTE Observations of IE 1547.0-5408 Prior to 2009 January 22 



Obs. Obs. Date MJD Num. PCUs Obs. Usage c Num. Profile 

Num. ID of On a Lcngth b of Codc c 

Bursts (ks) TOAs d 



1 


D93017-10-01-00* 


10/4/2008 


54743.36 


10-20 


2,3 


3.80 


t,p,f 


2 


A 


2 


D93017-10-01-02 


10/5/2008 


54744.41 





0,2 


0.59 


t 


1 


— 


3 


D93017-10-01-01 


10/5/2008 


54744.54 





2 


0.58 


t 


1 


_ 


4 


D93017-10-01-03 


10/6/2008 


54745.07 


0-10 


0,1,2,4 


1.40 


t,p,f 


2 


B 


5 


D93017-10-01-04 


10/6/2008 


54745.98 





1,2 


1.77 


t,p,f 


2 


C 


6 


D93017-10-01-08 


10/7/2008 


54746.31 


0-10 


1,2 


3.02 


t,p,fa 


2 


c 


7 


D93017-10-01-09 


10/7/2008 


54746.37 





2 


3.51 


t,p,fa 


2 


c 


8 


D93017-10-01-05 


10/7/2008 


54746.83 





2 


1.12 


t 


1 


— 


9 


D93017-10-01-ll g 


10/8/2008 


54747.03 





0,2 


2.19 


t,p,f 


2 


c 


10 


D93017-10-01-06 


10/8/2008 


54747.64 





2 


1.01 


t 


1 


— 


11 


D93017-10-01-10 


10/9/2008 


54748.92 





2 


1.65 


t,p,f 


1 


c 


12 


D93017-10-02-00 


10/10/2008 


54749.60 


0-10 


2,3,4 


1.12 


t,p,fa 


1 


c 


13 


D93017-10-02-01 


10/10/2008 


54749.67 





2,3p,4p 


1.63 


t,p,fa 


1 


c 


14 


D93017-10-02-02 


10/11/2008 


54750.10 





2 


2.73 


t 


1 




15 


D93017-10-02-06 


10/12/2008 


54751.21 





2,4 


0.93 


t,p,f 


1 


( 


16 


D93017-10-02-03 h 


10/13/2008 


54752.13 





0,2 


2.38 


t,p,f 


2 


D 


17 


D93017-10-02-04 


10/14/2008 


54753.37 





1,2 


3.46 


t,p,f 


2 


E 


18 


D93017- 10-02-05 


10/16/2008 


54755.33 





1,2 


3.49 


t,P,f 


2 


F 


19 


D93017-10-03-00 


10/17/2008 


54756.31 





2,4 


3.46 


t,p,f 


2 


G 


20 


D93017-10-03-01 


10/18/2008 


54757.58 





0,2 


1.10 


t,p,fa 


1 


H 


21 


D93017-10-03-02 


10/18/2008 


54757.62 





0,2 


1.21 


t,p,fa 


1 


H 


22 


D93017-10-04-00 


10/20/2008 


54759.00 





0p,lp,2,4p 


7.61 


t,p,f 


2 


H 


23 


D93017-10-04-01 


10/21/2008 


54760.96 





2,4 


3.49 


t,p,fa 


2 


H 


24 


D93017-10-04-02 


10/22/2008 


54761.37 





2 


1.53 


t,p,fa 


1 


H 


25 


D93017-10-04-03 


10/23/2008 


54762.20 


0-10 


2,3 


3.52 


t,p,fa 


2 


H 


26 


D93017-10-04-04 


10/23/2008 


54762.48 





2,4 


1.24 


t,p,fa 


1 


H 


27 


D93017-10-05-00 


10/24/2008 


54763.18 





2,4 


3.50 


t,p,f 


2 


H 


28 


D93017-10-05-02 


10/25/2008 


57764.84 





2,3 


1.05 


t,p,f 


1 




29 


D93017-10-05-03 


10/27/2008 


54766.06 





2,4 


2.06 


t,p,f 


2 




30 


D93017-10-05-01 


10/29/2008 


54768.80 





2 


3.04 


t,p,f 


2 




31 


D93017-10-06-00 


11/1/2008 


54771.81 





2,4 


3.92 


t,p,f 


2 




32 


D93017-10-06-01 


11/5/2008 


54775.02 





1,2 


3.69 


t,p,f 


2 




33 


D93017-10-06-03 


11/6/2008 


54776.98 





2,4 


1.66 


t,p,f 


2 




34 


D93017-10-07-00 


11/9/2008 


54779.99 





2 


1.65 


t,f 


2 




35 


D93017-10-07-01 


11/13/2008 


54783.07 





2,4 


1.48 


t,p,f 


2 


J 


36 


D93017-10-08-00 


11/16/2008 


54786.24 





2,3 


2.08 


t,p,f 


2 


.1 


37 


D93017-10-08-01 


11/20/2008 


54790.08 





2,4 


3.21 


t,p,f 


2 


.1 


38 


D93017-10-09-00 


11/23/2008 


54793.09 





2,3p 


6.26 


t,P,f 


2 


J 


39 


D93017-10-10-00 


12/3/2008 


54803.00 





1,2 


1.95 


t,f 


2 




40 


D93017-10-10-01 


12/3/2008 


54803.10 





2,3 


2.92 


t,p,f 


2 


J 


41 


D93017-10-10-02 


12/3/2008 


54803.20 





1,2 


1.92 


t 


2 




42 


D93017-10-11-00 


12/9/2008 


54809.90 





2,4p 


6.97 


t,p,f 


2 


J 


43 


D93017-10-12-00 


12/19/2008 


54819.60 





0p,2 


7.50 


t,p,f 


2 


K 


44 


D93017-10-12-01 


12/25/2008 


54825.70 





0p,lp,2 


8.34 


t,p,f 


2 


K 
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Table 1 — Continued 



Obs. 


Obs. 


Date 


MJD 


Num. 


PCUs 


Obs. 


Usage c 


Num. 


Profile 


Num. 


ID 






of 


On a 


Lcngth b 




of 


Code c 










Bursts 




(ks) 




TOAs d 




45 


D93017-10-13-00 


1/1/2009 


54832.60 





0,2 


6.15 


t,p,f 


2 


K 


46 


D93017-10-14-00 


1/8/2009 


54839.16 





2 


1.41 


t 


1 




47 


D93017-10-14-01 


1/8/2009 


58939.20 





2 


3.99 


t,p,f 


1 


K 


48 


D93017-10-15-01 


1/13/2009 


54844.97 





2,4 


2.73 


t,p,f 


2 


K 


19 


D93017-10-15-00 


1/14/2009 


54845.04 





2,3 


3.32 


t,p,f 


2 


K 



a If a PCU is on during one part of the observation only, its number is followed by a "p." Data from all PCUs were used in 
the timing analysis. Data from PCUs 2, 3, and 4 only were used in the pulsed flux analysis. 

b Total time on source. 

c "t" refers to this observation being used in the timing analysis, "p" refers to this observation being used in the pulse profile 
analysis, "f" refers to this observation being used in the pulse flux analysis, "fa" indicates two consecutive observations where 
the pulsed flux was averaged. 

d The number of independent TOAs extracted from this observation. 

c The letters in this column refer to the the time intervals indicated in the bottom panel of Figure [5] 

'There was one RXTE observation taken a few hours prior to D93017-10-01-00, but it was in a data mode different from all 
other obtained observations. 

g D93017-10-01-ll was originally named D93017-10-01-07. 

h The pulsed flux in this observation is significantly higher than that in the previous observations. 
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Table 2. RXTE Observations of IE 1547.0-5408 from 2009 January 22 to December 26 



Obs. Obs. Date MJD Num. PCUs Obs. Usage c Num. Profile 

Num. ID of On a Length b of Code 

Bursts (ks) TOAs d 



50 


D93017-10-17-00 


1/22/2009 


54853.86 


>20 


2 


6.66 


t,p,f 


2 


L 


51 


D93017-10-16-01 


1/23/2009 


54854.78 


10-20 


2,3 


2.64 


t,p,fa 


1 


L 


52 


D93017-10-16-00 


1/23/2009 


54854.84 


>20 


2,4 


3.35 


t,p,fa 


2 


L 


53 


D93017-10-16-05 


1/24/2009 


54855.82 


>20 


2 


10.29 


t,p,f 


2 


L 


54 


D93017-10-16-02* 


1/25/2009 


54856.81 


>20 


2 


5.97 


t,p,f 


2 


L 


55 


D94017-09-01-02 


1/26/2009 


54857.15 


>20 


2,4 


1.81 


t,p 


2 


L 


56 


D94017-09-01-01 


1/28/2009 


54859.88 


>20 


lp,2,3p 


10.32 


t,p,f 


2 


L 


57 


D94017-09-01-03 


1/29/2009 


54860.74 


>20 


0p,l,2,4p 


5.58 


t,p,f 


2 


L 


58 


D94017-09-02-01 


1/31/2009 


54862.56 


10-20 


0,2 


1.76 


t,p,fa 


1 


L 


59 


D94017-09-02-00 


1/31/2009 


54862.69 


10-20 


0p,lp,2,4 


17.46 


t,p,fa 


2 


L 


60 


D94017-09-02-02 


2/1/2009 


54863.75 


10-20 


lp,2 


5.91 


t,p,f 


2 


L 


61 


D94017-09-02-03 


2/1/2009 


54863.87 


>20 


2,4 


3.27 


t,p,fa 


2 


L 


62 


D94017-09-02-04 


2/1/2009 


54863.93 


>20 


0,2 


3.26 


t,p,fa 


2 


L 


63 


D94017-09-03-00 


2/6/2009 


54868.72 


10-20 


lp,2,3p 


5.81 


t,p,fa 


2 


M 


64 


D94017-09-03-01 


2/6/2009 


54868.83 


0-10 


2,4 


2.31 


t,p,fa 


1 


M 


65 


D94017-09-04-00 


2/13/2009 


54875.58 


10-20 


2,3p,4p 


7.61 


t,p,f 


2 


M 


66 


D94017-09-04-01 


2/18/2009 


54880.7 


0-10 


lp,2 


5.12 


t,p,f 


2 


M 


67 


D94017-09-05-00 


2/28/2009 


54890.1 


0-10 


1,2 


4.12 


t,p,fa 


2 


M 


68 


D94017-09-05-01 


2/28/2009 


54890.17 


0-10 


2,3 


1.58 


t,p,fa 


1 


M 


69 


D94017-09-06-00 


3/6/2009 


54896.84 





0p,lp,2 


7.24 


t,p,f 


2 


M 


70 


D94017-09-07-02 


3/16/2009 


54906.45 





1,2 


1.66 


t,p,fa 


2 


M 


71 


D94017-09-07-01 


3/16/2009 


54906.52 


0-10 


2,4 


2.01 


t,p,fa 


2 


M 


72 


D94017-09-07-00 


3/16/2009 


54906.59 





1,2 


2.37 


t,p,f 


2 


M 


73 


D94017-09-08-00 


3/21/2009 


54911.76 


10-20 


0p,lp,2 


5.87 


t,p,f 


2 


M 


74 


D94017-09-09-01 


3/30/2009 


54920.53 





2 


2.46 


t,p,f 


2 


M 


75 


D94017-09-09-00 s 


3/30/2009 


54920.59 


>20 


1,2 


3.06 


t,p,f 


2 


M 


76 


D94017-09-10-00 


4/5/2009 


54926.54 


0-10 


0p,lp,2 


5.84 


t,p,f 


2 


M 


77 


D94017-09-11-00 


4/14/2009 


54935.25 





2,4 


5.98 


t,p,f 


2 


M 


78 


D94017-09-12-00 


4/19/2009 


54940.87 


0-10 


lp,2,4p 


6.54 


t,p,f 


2 


N 


79 


D94017-09-13-00 


4/25/2009 


54946.63 





2,3 


2.98 


t,p,fa 


2 


N 


80 


D94017-09-13-01 


4/25/2009 


54946.7 


10-20 


2,4 


2.97 


t,p,fa 


2 


N 


81 


D94017-09-14-00 


5/3/2009 


54954.8 





2,4p 


6.98 


t,p,f 


2 


N 


82 


D94017-09-15-00 


5/11/2009 


54962.52 





2,4p 


5.95 


t,p,f 


2 


N 


83 


D94017-09-16-00 


5/18/2009 


54969.26 





2,3 


2.91 


t,p,fa 


2 


N 


84 


D94017-09-16-01 


5/18/2009 


54969.32 





2,4 


3.03 


t,p,fa 


2 


N 


85 


D94017-09-17-00 


5/26/2009 


54977.44 





2,3p,4p 


6.12 


t,p,f 


2 


N 


86 


D94017-09-18-00 


6/4/2009 


54986.19 





2,4p 


5.95 


t,p,f 


2 


N 


87 


D94017-09-19-02 


6/9/2009 


54991.83 





2,3 


1.48 


t,p,f 


2 


N 


88 


D94017-09-19-00 


6/10/2009 


54992.02 





2,3 


2.97 


t,p,fa 


2 


N 


89 


D94017-09-19-01 


6/10/2009 


54992.09 





2 


1.47 


t,p,fa 


1 


N 


90 


D94017-09-20-00 


6/17/2009 


54999.54 





2,3p,4p 


7.14 


t,p,f 


2 


N 


91 


D94017-09-21-00 


6/22/2009 


55004.6 





2,4 


5.6 


t,p,f 


2 


N 


92 


D94427-01-01-00 


6/30/2009 


55012.43 





2 


5.93 


t,p,f 


2 


N 


93 


D94427-01-02-00 


7/6/2009 


55018.26 





2,3p,4p 


5.89 


t,p,f 


2 
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Table 2 — Continued 



Obs. 


Obs. 


Date 


MJD 


Num. 


PCUs 


Obs. 


Usage c 


Num. 


Profile 


Num. 


ID 






of 
Bursts 


On a 


Length b 
(ks) 




of 
TOAs d 


Code 6 


94 


D94427-01-03-00 


7/14/2009 


55026.02 





2,3 


4.12 


t,p,f 


2 


O 


95 


D94427-01-04-00 


7/23/2009 


55035.66 


0-10 


2,3p,4p 


6.18 


t,p,f 


2 


O 


96 


D94427-01-05-00 


7/31/2009 


55043.44 





0p,lp,2 


7.32 


t,p,f 


2 


O 


97 


D94427-01-06-00 


8/11/2009 


55054.23 





0p,2,3p 


5.52 


t,p,f 


2 


O 


98 


D94427-01-07-00 


8/20/2009 


55063.2 





1,2 


1.89 


t,p,fa 


2 


O 


99 


D94427-01-07-01 


8/20/2009 


55063.25 


0-10 


2,3 


3.21 


t,p,fa 


2 


o 


100 


D94427-01-08-00 


8/28/20009 


55071.1 





0p,2 


3.35 


t,p,fa 


2 


o 


101 


D94427-01-08-01 


8/28/2009 


55071.16 





2 


3.63 


t,p,fa 


2 


o 


102 


D94427-01-09-00 


9/6/2009 


55080.99 





2 


7.72 


t,p,f 


2 


o 


103 


D94427-01-13-00 


10/8/2009 


55112.11 





0p,lp,2 


5.94 


t,p,f 


2 


p 


104 


D94427-01-14-00 


10/16/2009 


55120.94 





2,3 


3.63 


t,p,f 


2 


p 


105 


D94427-01-14-01 


10/17/2009 


55121.02 





2,4 


2.31 


t,p,f 


2 


p 


106 


D94427-01-15-00 


11/3/2009 


55138.14 





0p,lp,2 


5.95 


t,p,f 


2 


p 


107 


D94427-01-16-00 


11/11/2009 


55146.14 





0p,lp,2 


5.95 


t,p,f 


2 


p 


108 


D94427-01-10-00 


11/20/2009 


55155.22 





2,3 


5.93 


t,p,f 


2 


p 


109 


D94427-01-17-00 


11/27/2009 


55162.88 





2 


5.95 


t,p,f 


2 


p 


110 


D94427-01-11-00 


12/5/2009 


55170.92 





2 


6.24 


t,p,f 


2 


p 


111 


D94427-01-18-01 


12/12/2009 


55177.92 





2,4 


2.43 


t,p,fa 


2 


p 


112 


D94427-01-18-00 


12/13/2009 


55178.04 





2 


2.98 


t,p,fa 


1 


p 


113 


D94427-01-12-00 


12/19/2009 


55184.93 





lp,2,4p 


5.18 


t,p,f 


2 


p 


111 


D94427-01-19-00 


12/26/2009 


55191.12 





2,3p,4p 


6.38 


t,p,f 


2 


p 



a If a PCU is on during one part of the observation only, its number is followed by a "p." Data from all PCUs were used in 
the timing analysis. Data from PCUs 2, 3, and 4 only were used in the pulsed flux analysis. 

b Total time on source. 

c "t" refers to this observation being used in the timing analysis, "p" refers to this observation being used in the pulse profile 
analysis, "f" refers to this observation being used in the pulse flux analysis, "fa" indicates two consecutive observations where 
the pulsed flux was averaged. 

d The number of independent TOAs extracted from this observation. 

c Thc letters in this column refer to the time intervals indicated in the bottom panel of Figure [5] 
'ID D93017-10-16-02 is also called D94017-09-01-00. 

s In observation D94017-09-09-00, the count rate is unusually high, and one burst has a particularly long tail. 
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Table 3. Swift Observations of IE 1547.0-5408 from 
2008 October to 2009 December Used in this Paper 



Obs. ID 


Date 




MJD 


Obs. Length a 
(ks) 


Phase- 
Connected 10 


00330353001 


2008- 


■10- 


-03 


54742. 


.80 


14.9 


yes 


00330353002 


2008- 


■10- 


■04 


54743. 


79 


4.9 


yes 


00330353004 


2008- 


■10- 


-05 


54745. 


.14 


10.5 


yes 


00330353005 


2008- 


■10- 


-07 


54746. 


51 


7.7 


yes 


00330353006 


2008- 


■10- 


-08 


54747. 


11 


4.5 


yes 


00330353007 


2008- 


•10- 


-09 


54748. 


39 


3.6 


yes 


00330353008 


2008- 


■10- 


-10 


54749. 


.46 


3.8 


yes 


00330353010 


2008- 


■10- 


-12 


54751. 


.53 


3.9 


yes 


00330353011 


2008- 


■10- 


•13 


54752. 


.40 


3.4 


yes 


00330353012 


2008- 


■10- 


-16 


54755. 


11 


1.0 


yes 


00330353013 


2008- 


•10- 


•17 


54757. 


.38 


5.0 


yes 


00330353014 


2008- 


10- 


-20 


54759. 


.80 


3.9 


yes 


00330353015 


2008- 


■10- 


-22 


54761. 


.43 


3.8 


yes 


00330353016 


2008- 


■10- 


-24 


54763. 


21 


3.5 


VPS 


00090007024 


2009- 


-01- 


•04 


54835. 


.13 


3.3 


no 


00090007025 


2009- 


■01- 


-12 


54844. 


.13 


4.2 


no 


00340573000 


2009- 


■01- 


-22 


54853. 


20 


7.6 


yes 


00340573001 


2009- 


-01- 


-22 


54853. 


55 


9.4 


yes 


00090007026 


2009- 


■01- 


-23 


54855. 


05 


8.2 


yes 


00090007027 


2009- 


■01- 


-25 


54856. 


19 


3.2 


yes 


00090007028 


2009- 


-01- 


-26 


54857. 


.25 


3.5 


yes 


00090007035 


2009- 


-02- 


•02 


54864. 


72 


4.1 


yes 


00090956040 


2009- 


■02- 


•05 


54867. 


.57 


6.1 


yes 


00090007036 


2009- 


-02- 


•12 


54874. 


32 


6.0 


yes 


00090007037 


2009- 


-02- 


-22 


54884. 


.63 


4.6 


no 


00090007038 


2009- 


■03- 


•04 


54894. 


.61 


3.9 


no 


00090007039 


2009- 


03- 


•13 


54904. 


.45 


4.0 


no 


00090007040 


2009- 


03- 


•24 


54914. 


.90 


4.1 


no 


00030956043 


2009- 


■04- 


•29 


54950. 


.47 


1.6 


no 


00030956045 


2009- 


■05- 


•27 


54978. 


.52 


2.1 


no 


00030956053 


2009- 


-09- 


•16 


55090. 


52 


1.4 


no 


00030956054 


2009- 


■09- 


•30 


55104. 


.63 


3.3 


no 


00030956055 


2009- 


■10- 


•14 


55118. 


09 


1.9 


no 



a XRT time on source. 

b "yes" means that this observation belongs to a period of time covered by a 
unique phase-connected timing solution reported in this paper. See Section [3] 
in the text. 
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Table 4. Spin Parameters for IE 1547.0— 5408 a for the 
Phase-Connected Segments 



Parameter 


Post-2008 outburst 
ephemeris spanning 
29 days 


Post-2009 outburst 
ephemeris spanning 
23 days 


MJD range 


54743.4-54771.8 


54853.9-54875.7 


No. TOAs b 


46 RXTE + 14 Swift 


31 RXTE + 8 Swift 


v (Hz) 


0.48277893(4) 


0.48259625(3) 


v (1(T 12 Hz s- 1 ) 


-6.19(8) 


-5.21(4) 


v (10- 18 Hz s- 2 ) 


-6.69(7) 


< l.l c 


Epoch (MJD) 


54743.000 


54854.000 


RMS residual 


3% 


4% 



a Numbers in parentheses are TEMPO-reported la uncertainties. 

b The number of RXTE TOAs is larger than the number of 
RXTE observations because the longest observations yielded mul- 
tiple TOAs. 



c 3<7 upper limit. 



Year 

2009 2010 



J i I i u i i i I i i i i i i I I i i i I 

54800 54900 55000 55100 55200 

Modified Julian Day 



Fig. 1. — Epochs of the observations of IE 1547.0—5408 used in this paper. The first row 
of vertical bars marks the location of RXTE observations. The second row of vertical bars 
marks the location of Swift observations. The two solid vertical lines indicate the onset of 
the two outbursts. 
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Year 



10 — 



0.4832 
0.4830 
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Fig. 2. — Long-term evolution of IE 1547.0—5408 properties, (a) Fre quency as a function of 
time. Green square: frequency from radio data ( ICamilo et al.l 120071 ). Green li ne: evolution 



of the frequency u as a quadratic polynomial where v Q and v Q are o btained fromlCamilo et al. 
( 20071 ). and v is from the best-fit slope to Panel 2 of Figure 6 in ICamilo et al. (j2008|). Red 
solid lines: frequency as a function of time from phase-coherent timing. Red dotted lines: 
possible but not unique phase- coherent solutions. Dark blue points: individual frequency 
measurements, each obtained from a combination of closely spaced observations. Light blue 
points: individual frequency measurements, each obtained from a long observation, (b) 
Frequency de r ivativ e as a function of time. Green points: z> measurements taken from 
Camilo et al.l (120081 ). Red solid lines: z> as a function of time from phase- coherent timing. 
Red dotted lines: possible but not unique phase-coherent solutions, (c) Pulsed flux in the 
2—10 keV band, (d) Count rate in the 2—10 keV band from the full RXTE field of view. 
All panels: the two solid vertical lines mark the onset of the two outbursts. 
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-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 

v( xl0~ u s~ 2 ) 



Fig. 3. — Results of the search for a phase-coherent timing solution consisting of a frequency 
v and of a single frequency derivative v for a set of TOAs covering a time period of 11 days, 
and starting 21 days after the onset of the first outburst (MJD 54764.846 to 54775.050). 
The x-axis shows trial v values. The y-axis shows trial v values. For every pair of v and z>, 
the shading corresponds to the reduced \ 2 obtained when the pair was used as set of trial 
phase-coherent timing parameters by TEMPO. White squares correspond to a reduced \ 2 
of 6 and above. Black squares correspond to a reduced x 2 of 0. The center of each of the 
uniformly-shaded islands corresponds to the best-fit v and z> found by TEMPO for each point 
on that island, and the size of each island corresponds to the uncertainty on the obtained v 
and v. The presence of several black islands clearly indicates the presence of several plausible 
timing solutions differing from each other by a small number of phase counts. 
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0.4B22 - - 

54800 54900 55000 55100 55200 

Modified Julian Day 

Fig. 4. — Frequency as a function of time for IE 1547.0—5408. Shaded black columns: each 
of the black columns corresponds to one of the overlapping sets of TOAs. The x-coordinate 
of the column corresponds to the time in the middle of that set of TOAs. Every column 
was obtained by collapsing vertically a rectangle similar to that shown in Figure [3j The 
shading of each square in every shaded column appearing in this Figure corresponds to 
the darkest shading in the corresponding row in a v and i> map similar to that shown in 
Figure [31 and represents the best-fit reduced x 2 - The presence of several black squares in 
every column indicates the presence of several plausible phase-coherent timing solutions at 
each epoch. The black horizontal lines represent the limits of the searched frequency values. 
The red solid lines represent frequency as a function of time from phase- coherent timing. 
The red dotted lines represent possible but not unique phase-coherent solutions. The dark 
blue points represent individual frequency measurements, each obtained from a set of closely 
spaced observations. The light blue points represent individual frequency measurements, 
each obtained from a long observation. The dashed line is an eyeball fit to all frequency 
measurements from 2008 and 2009. The two vertical black lines denote the onset of the two 
outbursts. 
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Fig. 5.— Normalized 2-6.5 keV pulse profiles of IE 1547.0-5408 from 2008 October to 
2010 December. Two cycles per plot are shown for clarity. The letter shown in the top- left 
corner of each plot refers to the time segments marked by arrows in the bottom plot, where 
the 2—10 keV RMS pulsed flux is shown for reference. The profiles marked with the letters 
A, B, and D to G are obtained from individual observations. The first and third dotted lines 
in the bottom plot correspond to the onsets of the two outbursts. 
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Fig. 6.— Normalized 2-4 keV, 4-10 keV, and 2-10 keV pulse profiles of IE 1547.0-5408 
for two selected time segments showing the dependence of the location of the peak on energy. 
Two cycles per plot are shown for clarity. 
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Fig. 7. — Pulsed flux time series and pulsed hardness ratio of IE 1547.0—5408. Top: evolu- 
tion of the RMS pulsed flux of IE 1547.0-5408 in 2-4 keV, 4-10 keV, and 2-10 keV. Note 
the latter time series was shown by Ng et al. (2011) and Scholz & Kaspi (2011). Bottom: 
The hardness ratio as a function of time, computed from the ratio of pulsed fluxes in the 
energy range 4—20 keV to 2—4 keV. All panels: The first dotted line marks the onset of the 
first outburst. The second dotted line marks the location of an enhancement in the pulsed 
flux. The third line marks the onset of the second outburst. The last two dotted lines bracket 
the time interval inside which the hardness ratio rose. 
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Fig. 8. — Comparison of the timing and pulsed flux properties of IE 1048.1—5937 (left 
column) and IE 1547.0—5408 (right column), (a) Frequency as a function of time for 
IE 1048.1— 593 7 with a linear trend subtracted (for details, see Figures 3 and 15 of 



Dib et al.l 120091 ). (b) Frequency as a function of time for IE 1547.0—5408 with a linear 
trend subtracted (see Figure 2 for details), (c) Frequency derivative as a function of time 
for IE 1048.1—5937. (d) Frequency derivative as a function of time for IE 1547.0—5408. (e) 
Pulsed flux as a function of time in the 2—10 keV band for IE 1048.1—5937. (f) Pulsed flux 
as a function of time in the 2—10 keV band for IE 1547.0—5408. 
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